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Abstract 

Gaussian sum-rules relate a QCD prediction to a two-parameter Gaussian-weighted 
integral of a hadronic spectral function, providing a clear conceptual connection to quark- 
hadron duality. In contrast to Laplace sum-rules, the Gaussian sum-rules exhibit enhanced 
sensitivity to excited states of the hadronic spectral function. The formulation of Gaussian 
sum-rules and associated analysis techniques for extracting hadronic properties from the 
sum-rules are reviewed and applied to scalar gluonium. With the inclusion of instanton 
effects, the Gaussian sum-rule analysis results in a consistent scenario where the gluonic 
resonance strength is spread over a hroad energy range below 1.6 GeV, and indicates the 
presence of gluonium content in more than one hadronic state. 
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1 Introduction 



The hadronic spectrum has too many scalar states above 1 GeV for a qq nonet, as would 
be anticipated if gluonium states exist in the 1-2 GeV region Determining how 
this gluonium content is distributed among these scalar-isoscalar resonances is thus an 
important issue. In particular, the possibility that the observed hadronic states are mixtures 
of gluonium and quark mesons must be explored. 

Gaussian sum-rules are sensitive to the hadronic spectral function over a wide energy 
range, and are thus well-suited to studying the distribution of gluonium states. The 
simplest Gaussian sum-rule (GSR) 

to ^ ^ 

relates a QCD prediction on the left-hand side of (|lD to the hadronic spectral function 
pit) (with physical threshold to) smeared over the energy range s — 2-\/r < t < s + 2-\/r, 
representing an energy interval for quark-hadron duality. An interesting aspect of the GSR 
is that the duality interval is constrained by QCD. A lower bound on this duality scale r 
necessarily exists because the QCD prediction has renormalization-group properties that 
reference running quantities the the energy scale z/^ = ^/t [0, |^. Thus it is not possible 
to achieve the formal r — >^ limit where complete knowledge of the spectral function 
could be obtained via 

\imG{s,r) = -p{s) , s > to . (2) 

T^O TT 

However, there is no theoretical constraint on the quantity s representing the peak of 
the Gaussian kernel appearing in (|l|). Thus the s dependence of the QCD prediction 
G{s,t) probes the behaviour of the smeared spectral function, reproducing the essential 
features of the spectral function. In particular, as s passes through t values corresponding 
to resonance peaks, the Gaussian kernel in (|l]) reaches its maximum value, implying that 
Gaussian sum-rules weight excited and ground states equally. This behaviour should be 
contrasted with Laplace sum-rules 



oo 




to 



where excited states are damped by the exponential decay of the Laplace kernel. Thus, 
in comparison with Laplace sum-rules the Gaussian sum-rule (GSR) has an enhanced 
sensitivity to excited states of the spectral function. 

In this paper, the original formulation of Gaussian sum-rules and analysis techniques 
for extracting spectral function hadronic features [|], ^ will be reviewed. These techniques 
will then be applied to scalar gluonium, where instanton contributions are known to be 
crucial for a consistent Laplace sum-rule analysis ||. Results of the GSR analysis 
indicate that the gluonium spectral strength is distributed across a broad energy range 
below 1.6 GeV [§. 
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2 Foundations of Gaussian Sum-Rules 



The general formulation of GSRs will be reviewed in the context of scalar gluonium probed 
by the following correlation function. 



n(Q2) = i J d^a;e'^-^(0|T{ J(x), J(0)}|O) , = -q' 



(4) 



^2 



11 1 ^51 19 

/5o = ^-gn; , A = y-^^/ , ••■ (7) 

The current J{x) is renormalization-group invariant in the chiral limit of Uf massless 
quarks as needed to probe physical (renormalization-group invariant) hadronic states. 

From the asymptotic form and assumed analytic properties of (§) follows a dispersion 
relation with three subtraction constants 

n(Q2) _ n(o) = g2n'(o) + ^Q^n"(o) -91 T p^A^^t , q''>o (8) 



2^ ' ' 71 4 t^{t + Q 

where p(t) is the hadronic spectral function with physical threshold to. and the subtraction 
constant 11(0) has been included on the side of the equation containing the QCD prediction 
because it is determined by the low-energy theorem [0] 

n(0)^ lim n(Q2) = |^(J) . (9) 

The undetermined subtraction constants n'(0) and n"(0) and field-theoretical divergences 
in n(Q^) proportional to Q'^ are eliminated in an integer-weighted family of Gaussian sum- 
rulesQ 

G,ii,r) . { '-^ + '^)^"(-^ - - <^ - '^>^"(-^ + } (10) 

where k = —1, 0, 1, . . . and with the Borel transform B defined by 

(-A2)^ / d \^ 

A2/Af=4T 

Applying definition ( |10D to both sides of (|) annihilates the undetermined low-energy 
constants and the field theoretical divergence contained in Il(Q'^). Using the identity 



B 



(A 



2\n 



A2 + a 



^(-a)"exp('^ ) forn>0 , (12) 



"'^This definition is a natural generalization of that given in Q|. To recover the original Gaussian 
sum-rule, we simply let fc = in (lir 



3 



then leads to the following GSR family 

^2 



Gk(s,T) + Sk-1 



41 



: exp 



vrr 



4r 



n(o) 



exp 



to 



4r 



TT 



-p(t)dt (13) 



Calculation of the Borel transform is achieved through an identity relating ([Tl| ) to the 
inverse Laplace transform H 



i3[/(A^)] = -/:-n/(A^)] 



where, in our notation, 



1 /A2\ 



(14) 



(15) 



with a chosen such that all singularities of / lie to the left of a in the complex A^-plane. 
With a change of variables, the calculation of the GSR reduces to [p[ 



Gfc(s,r) 



1 1 



A/iTrr 27ri Jn+ra 



-w) exp 



-(s + w)^ 



4r 



n(w)dw 



where Fi and r2 are the parabolae depicted in Figure |l|. 




(16) 




Figure 1: Contour of integration V\ + r2 defining the Gaussian sum-rule. The wavy hue 
on the negative real axis denotes the branch cut of n(z). 



3 Normalized Gaussian Sum-Rules 

Studies of Gaussian sum-rules have traditionally focussed on their connection with finite- 
energy sum-rules as established through the diffusion equation 

dp Or ' ^ ' 
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In particular, the resonance(s) plus continuum model 

p(t)=p^^'^(t)+^^(t-So)ImnQC°(t) , 



(1^ 



when p^^^{t) is evolved through the diffusion equation (|17D, only reproduces the QCD 
prediction at large energies (r large) if the resonance and continuum contributions are 
balanced through the finite-energy sum-rules [|| 



so 



n = integer 



(19) 



to 



Within the resonance(s) plus continuum model ([T^ ), the continuum contribution to 
the GSRs is determined by QCD 



V4: 



TCT 



exp 



4r 



TT 



-ImnQC°(t)dt 



and is thus combined with Gk (■§, r) to give the total QCD contribution 

° (s, r, so) = Gk (5, r) - G^ {s, r, s,) , 
resulting in the final relation between the QCD and hadronic sides of the GSRs. 

At 



(20) 



(21) 



° {s, T, So) + 6k -1^= exp [ ) n(0) = H exp 
V47rr V 4r / Jt^, 



TT 



(22) 

Integrating both sides of (|2^ ) reveals that the normalization of the GSRs is related to 
the finite-energy sum-rules 



oo oo 

j G?^°(s,r,so)ds + 4_in(0) = ^ J tV'Wdt 



(23) 



to 



Thus the diffusion equation analysis relates the normalization of the GSR to the 
finite-energy sum-rules. Information independent of this relation is extracted from the 
normalized GSRs 



(24) 
(25) 



_ ^r(g,r,^o) + 4-i^exp(^ )n(0) 
' M2J°(r, .o) + 4-in(0) 

oo 

Mk,n{T,So)= j s^Gk{s,T,so)<ls , n = 0, 1, 2, . 



which are related to the hadronic spectral function via 



1 TOO ,h 

L ^ exp 



QCD, 



s,r, So) 



/47rr J to 



-js-tf 

4r 



(26) 
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4 Gaussian Sum-Rule Analysis Techniques 



In the single narrow resonance model, p^^'^{t) takes the form 



(27) 



where m and / are respectively the resonance mass and coupling. With such an ansatz 



the normalized Gaussian sum-rule ( P6D becomes 

1 



: exp 



vrr 



At 



(28) 



Deviations from the narrow-width limit are proportional to m^F^/r, so this narrow-width 
model may actually be a good numerical approximation. Phenomenological analysis of the 
single narrow resonance model proceeds from the observation that (^) has a maximum 
value (peak) at s 



m 



independent of the value of r. The value of sq is then optimized 
by minimizing the r dependence of the s peak position of the QCD prediction, and the 
resulting r-averaged s peak position leads to a prediction of the resonance mass 0. 

The p meson illustrates this single-resonance analysis technique and demonstrates 
that GSRs can be used to predict resonance properties. The correlation function of the 
vector- isovector correlation function results in the following {k = 0) GSR 



G 



QCD 



S,T, Sq 



1 + 



a 



iVry\ 



TT 



erf 



erf 



So 



2v^ 



+ 



327r2r 
1 



exp - — 
TIT \ At 



GAn'^Ty/nr 
s 



-1 + 



2r 



2V^ 



%2 \ 



-i7 ) (cm) 



where 



{cm) = I {aC) - 8n'm {qq) 



(29) 
(30) 



and SU{2) symmetry along with the vacuum saturation hypothesis have been employed. 
For brevity, we refer to the literature for the expressions for the dimension eight 
condensates, and simply use ( ^Qf) to establish a convention consistent with [|^. Note that 
the running coupling is referenced to the energy scale ^/r, a point which will be discussed 
in the next Section. 

The non-perturbative QCD condensate contributions in (|2^) are exponentially sup- 
pressed for large s. Since s represents the location of the Gaussian peak on the phe- 
nomenological side of the sum-rule, the non-perturbative corrections are most important 
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in the low-energy region, as anticipated by the role of QCD condensates in relation to the 
vacuum properties of QCD. This explicit low-energy role of the QCD condensates clearly 
exhibited for the Gaussian sum-rules is obscured in the Laplace sum-rules. 
The QCD inputs used for the p meson analysis analysis are 



A(^) = 300 MeV 

G^) = (0.045 ± 0.014) GeV^ , 2m {qq) 
8Q6 

{ClOD = -fvs^n' (1.8 X lO-^GeV^) , U 
(Cg^O^) = (0.40 ±0.16) GeV^ , 



r2 2 
1.5 ±0.5 



(31) 
(32) 

(33) 
(34) 



consistent with the condensate parameters in 
predicts the following values for rup and sq 



The criteria of r stability of s peak 



(0.75 ± 0.07) GeV , Sq = (1.2 ± 0.2) GeV^ 



(35) 



in excellent agreement with the known value of the p mass. Furthermore, the phenomeno- 
logical and QCD sides of the (normalized) Gaussian sum-rules shown in Figure ^ are in 
superb agreement even for very small values of r 0. This agreement is particularly 
impressive since there are no free parameters corresponding to the normalization of the 
curves in this analysis. 

In more complicated resonance models, the s peak position of the phenomenological 
model begins to develop r dependence which is well-described by O, 



. B C 
A + — + — 

T T'^ 



(36) 



Analysis of how the s peak "drifts" with r in comparison with the behaviour ( ^6l) then 
allows optimization of sq. After optimization of sq, the resonance model parameters are 
extracted from various moments of the GSRs. For example, a square pulseQ centred at 
t = rri^ with total width 2mr leads to the following [k = 0) normalized GSR 



iVo^''°(5,r,.o) 



4mr 



erf 



s — ni^ ± mV 



erf 



mV 



2v^ 



(37) 



Expansion of (|^) for small V demonstrates that deviations from the narrow width limit 
(^) scale as m^F^/r. The resonance parameters can then be determined by the following 



moment combinations [see (|25| )] of the right-hand side of (|37|) 



m 



0,0 



2r = 



Mr 



0,0 



0,0 



2-p2 



-m r 



(38) 
(39) 



^This model would describe a broad, structureless feature in the hadronic spectral function. 
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Figure 2: Comparison of tlie vector-current theoretical prediction for N^'~'^ {s,t, sq) 
with the single resonance phenomenological model. The r values used for the four pairs 
of curves, from top to bottom in the figure, are respectively r = 0.5 GeV^, r = 1.0 GeV^, 
r = 2.5 GeV^, and r = 4.0 GeV^. Note the almost complete overlap between the QCD 
prediction and phenomenological model. 

where it is understood that the QCD expressions at the optimized value of sq are used on 
the left-hand side.0 

Moments also provide a method for testing the accuracy of agreement between the 
QCD and phenomenological sides of the normalized GSR beyond a simple measure 
which could be extracted from plots such as Figure ||. For example, a combination of 
third-order moments representing the asymmetry of the s dependence about its average 
value results in Wi 



A, 



(3) 



0,0 



Mo.2\ /Mo. 



M 



0,0 



0,0 



0,1 



Mr 



0,0 



(40) 



and hence a deviation of the QCD value of this moment from its value of A'^^ = in the 
square pulse model indicates a shortcoming of the model spectral function in comparison 
with QCD. 



■^The residual moment combinations leading to the resonance parameters must be r independent, 
providing a consistency check on the analysis. Weak residual r dependence is averaged over the t 
range used to extract the optimized sq from the peak-drift analysis. 
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The procedure for studying an A^-parameter resonance model is easily generalized. The 
peak-drift analysis is used to optimize Sq, the lowest moments are used to determine 
the resonance model parameters, and the next-highest moment combination is employed 
as a test of the accuracy of the model's agreement with the QCD prediction. 

5 Scalar Gluonium Gaussian Sum-Rules 

The lowest two Gaussian sum-rules for scalar gluonium contain perturbative, condensate 
and instanton corrections, and are given by H] 



G^^^{S,T, So) 



1 



rcT Jo 



t dt exp 



At 



[ao - TT aa) + 2ai log — 



t 



+3a2log 



1 / -s 



v47rr 

i67r3 r , , , . n 

—== / dn{p)p / texp 
V47rr J Jo 

1287r2 



V4 



exp 



7rr 



4r 



dn(p) 



At 



J2 (pVi) Y2 (pVt) dt 



(41) 



Go'^^is,T, So) 



\l A-KT Jo 



so 



t^dtexp 



[s-tl 
At 



1 r° 

—==hi{J) I exp 
\/47rr Jo 



At 



(ao - vr 02) + 2ai log ^ 



+3a2 log^ ( 4j 



+ 



1 



A/47rr 



exp 



4r 



TTT 



so 

dn{p)p'^ I exp 




-[s 



At 



J2 [pVtj Y2 i^pVtj dt . 

(42) 



The perturbative coefficients in these expressions are 



ao 



' I j) 



659 a , /ax ' 

— - + 247.480 - 
6b TT Vtt/ 



ai 



ax 3 



TT 



9 a 
- + 65.781- 

4 TT 



, aa = -10.1250 (^^) 



(43) 
(44) 
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as obtained from the three-loop MS calculation of the correlation function []TD|]. As a 
result of renormalization group scaling of the GSRs 0, the coupling in the perturbative 
coefficients is implicitly the running coupling at the scale z/^ = i/r in the MS scheme 



TT 



log 



A2 



A . _9 



Pi 



P2 



3863 
"38? 



(45) 
(46) 



L 

with Ajfg ^ 300 MeV for three active flavours, consistent with current estimates of 

«(M.) 0. 

The condensate contributions in (0), ( |42D involve next-to- lea ding order [0 contribu- 
tionsQfrom the dimension four gluon condensate (J) and leading order [|T3|] contributions 
from gluonic condensates of dimension six and eight 



(a) 



(47) 
(48) 



At: — 

TT 



1 + 



175 a 

36 TT 



hi 



Co 



c/o = 8 



TT^- . (49) 
TT 



The remaining terms in the GSRs (PT|), represent instanton contributions obtained 
from single instanton and anti-instanton []TT]] (z.e., assuming that multi-instanton effects 
are negligible |]T^) contributions to the scalar gluonic correlator |T^, |T5|. The 

quantity p is the instanton radius, n{p) is the instanton density function, and J2 and Y2 
are Bessel functions in the notation of JTBH. 

The instanton contributions to the Gaussian sum-rules can be interpreted [|| as nat- 
urally partitioning into an instanton continuum portion devolving from 



-Imm""*(t) 

TT 



-IQ-K^ j dn(p) pV J2 {pVPj Y2 [pVi) 



(50) 



and a contribution which, like the 11(0) low-energy theorem (LET) term, appears only in 
the k = —1 sum-rule 



1287r' 



exp 



—s 

17 



dn(p) 



(51) 



This asymmetric role played by the instanton is crucial in obtaining a consistent analysis 
from these two sum-rules. In the absence of instanton contributions the LET tends to 
dominate the left-hand side of (P^, corresponding to a massless state in the single-narrow 
resonance model (^81). However, the higher-weighted sum-rules are independent of the 
LET, and in the absence of instantons lead to a much larger mass scale in sum-rule 
analyses. This discrepancy between the LET-sensitive and LET-insensitive sum-rules has 



"'The calculation of next-to-leading contributions in |12 have been extended non-trivially to 
Uf = 3 from Uf = 0, and the operator basis has been changed from (aG^) to (J). 
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been shown to be resolved when instanton contributions are included in the Laplace sum- 
rules 0. A similar qualitative behaviour emerges from the Gaussian sum-rules, since the 
LET term in (^) and the LET-like instanton contribution (0) have the same functional 
form [i.e., each is proportional to exp (— s^/(4r))] and occur with opposite sign in the 
left-hand side of Thus there exists a cancellation between these effects, which is 



easily verified as being significant in the instanton liquid model 



n(p) = nj^p - Pc) 



8.0 X 10-* GeV* , Pc 
along with a standard value for the gluon condensate [IB 



0.6 



(aG^) = (0.07 ±0.01) GeV^ 



(52) 
(53) 

(54) 



This qualitative argument is upheld by the detailed GSR analysis presented in [Q]. However, 
the k = —1 GSR analysis is more sensitive to QCD uncertainties, justifying a focus on 
the A; = GSR for the remainder of this paper. 



6 Gaussian Sum- Rule Analysis of Scalar Gluonium 

A single narrow resonance model analysis of the k = Gaussian sum-rule results in a mass 
scale of approximately 1.3 GeV, but leads to poor agreement between the phenomeno- 
logical model and QCD prediction as shown in Figure 0, indicating that the gluonium 
spectral function is poorly described by a single narrow resonance. Furthermore, in the 
narrow width model the second-order moment combination should satisfy 

- 2r = , (55) 

but Figure ^ shows a substantial deviation from this behaviour [Q. Since this moment 
combination is related to the width of the GSR, we conclude that the gluonium resonance 
strength must be distributed over a significant energy range. 

The clear failure of the single narrow resonance model, indicative of distributed reso- 
nance strength significant enough to be resolved by the GSRs, is a significant conclusion in 
its own right, but various distributed resonance strength models have also been analyzed 
0. In order of increasing number of parameters (and increasing complexity) they are 

1. Single non-zero width models (2 parameters: mass, width) 

2. Two narrow resonance model (3 parameters: two masses, relative resonance strength) 

3. Narrow resonance plus a non-zero width resonance models (4 parameters: two 
masses, one width, relative resonance strength) 
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0.10 




Figure 3: Comparison of the theoretical prediction for Nq'^'^ {s,t,so) with the single 
narrow resonance phenomenological mode The r values used for the three pairs of curves, 
from top to bottom in the figure, are respectively r = 2.0 GeV^, r = 3.0 GeV^, and 
r = 4.0 GeV^. Note the prominent disagreement between the QCD prediction and phe- 
nomenological model in the vicinity of the peaks. 



In the single non-zero width models, elaborations on the square pulse include a Gaus- 
sian resonance and a skewed Gaussian resonance models 



p{t) ~ exp 



m 



2^2 



212 



Pit) 



exp 



{t — m 



2r2 



(56) 



(57) 



which are analytically and numerically simpler to analyze than a Breit-Wigner shape. In the 
Gaussian resonance models, the quantity T can be related to an equivalent Breit-Wigner 
width by Tbw = v^21og2r/m, and the factor in the skewed Gaussian is chosen 
for consistency with (low-energy) two-pion decay rates [|^, The relevant moment 
combinations for the Gaussian model ( ^6)) are [|| 

Mo,i 



where 



+ TA 

Mo,o 

(j^ - 2r = - m^TA - T^A^ 

= (m^r - r^) A + sm^r^A^ + 2T^a^ 



A 



exp -2r2 



(58) 

(59) 
(60) 

(61) 
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T (GeV*) 



Figure 4: Plot of ctq for the theoretical prediction (dotted curve) compared with cTq = 2r 
for the single-resonance model (solid curve) for the k = sum-rule. 

The quantity A is small, so the Gaussian model cannot accommodate large values of the 
asymmetry A^q \ However, the skewed Gaussian naturally leads to a larger asymmetry as 
reflected in the following results for the relevant moment combinations 

m2(m^ + 3r2) 



Mr 



0,1 



Mr 



0,0 



2r 



(m4 + r2)2 



+ 0(A) 



+ 0(A) 



A, 



(3) 



(m4 + r2)3 



+ 0{A) 



(62) 
(63) 
(64) 



The results for the single non-zero resonance models are shown in Table 0, and indicate 



that the non-zero width models underestimate the QCD value of the asymmetry Aq^-* by 
at least an order of magnitude 0. This failure suggests that further phenomenological 
models which can generate a larger asymmetry are required. 





mass (GeV) 


width (GeV) 


(GeV^) 


square pulse 


1.30 ±0.17 


0.59 ±0.07 





unskewed gaussian 


1.30 ±0.17 


0.40 ±0.05 


0.000342 


skewed gaussian 


1.17±0.15 


0.49 ±0.06 


0.00943 


QCD 






-0.0825 



Table 1: The results of a /c = Gaussian sum-rules analysis of a variety of non-zero 
resonance width models. The quoted resonance parameters include uncertainties introduced 
by the QCD input parameters, except for Aq ' which is obtained from the central values. 
For the Gaussian resonance models, the given width is the equivalent Breit-Wigner width 
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A model containing two narrow resonances 



results in the following normalized Gaussian sum-rule {I 



TCT 



ri exp 



m 



2\2 



4r 



+ r2 exp 



m 



2\2 



4r 



where 



ri 



r2 



/I 



f? + /I 



ri+r2 = 1 



(65) 



(66) 



(67) 



/l + /2 Jl ~^ J2 

parameterize the relative strength of the two resonances. This model can accommodate 
a large asymmetry parameter through an asymmetric distribution of resonance strength, 
as indicated by the following moment combinations 



Mo,i 



M 



'-(z + ry) 



0,0 



A, 



(3) 



2r=-t/'(l 



So - 12r' - 12r [a', - 2r) = (l - r^) (l + Sr^) 



where 



2 I 2 

z = nil + ^^2 ? 



^0 = 



0,4 



Mo,3Mo,i , ^Mo,2 /^Mo,i 
4^rT — h 6 



0,0 



Mo,o Mo,o 



Mr 



0,0 



0,0 



ri - r2 



Mo,i 



M 



0,0 



(68) 

(69) 
(70) 
(71) 

(72) 
(73) 



The lowest three moment combinations (|6^)-(|70|) are used to determine the three reso- 
nance parameters, and the final fourth-order residual combination ( [TTD is used as a test of 
the agreement between QCD and the phenomenological model. The resulting resonance 
parameters are 



mi 



= (0.98 ± 0.2) GeV, ms = (1.4 ± 0.2) GeV, n = 0.28 T 0.06, ra = 1 - n. (74) 

Results from the fourth-order moment measure of the accuracy between QCD and the 
two-narrow resonance model is given in Table 0, and indicate an approximately 50% 
discrepancy between QCD and the phenomenological model 0. 

As a final attempt to improve the agreement between the next-highest moment test 
of the agreement between QCD and the phenomenological model, consider a narrow 
resonance of mass m and relative strength r^, combined with a second resonance of 
mass M, width F and relative strength rj\/. Of course for a normalized GSR we must 
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So - 12r2 - 12r {a^ - 2r) 


QCD 


0.170 GeV« 


double resonance model 


0.074 GeV^ 



Table 2: The results of a A; = Gaussian sum-rules analysis of the double narrow resonance 
model using central values of the QCD parameters. 



have + rM = 1, so this defines a four-parameter model. For example, when a square 
pulse is used to describe the second resonance the resulting normalized GSR is 



N 



QCD 



(s, r, So) =rr, 



: exp 



vrr 



2\2 



(s — m ) 



+ Tm 



AMY 



erf 



4r 

s - m2 + Mr 



— erf 



M2 - MT 
2v^ 



(75) 



where rm + ^A/ = 1. Four moment combinations are then needed to define the resonance 
parameters, and a fifth-order moment combination serves as a measure of the agreement 
between QCD and the phenomenological model [H] 



M 



:z + r?/) + C(A) 



0,0 



- 2r = (1 _ ^2) ^ _^r2 _ ^) (1 _ ^) + (^) 
-^yV (1 - r') - ^T'y (1 - r') {z - y) + O (A) 



(3) 



(76) 

(77) 
(78) 



^0 - I2r' - 12r {a^ - 2r) =ly^ (l - r^) (l + 3r') + (1 + r) (l - r') {z - y) 

it) o 



(5) 



20rA 



(3) 



+ -rMl-r) (z-y)^ + 0(A) 



- (1 - r^) (1 + r^) - V (1 - r') (1 + r)^ (;. - y) 



(79) 



(80) 



where the fifth-order asymmetry moment is 
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M, 



k,0 



Ml 



k,0 



+ 4 



Ml 



k,0 



The resulting resonance parameters and the QCD prediction and phenomenological values 
of the fifth-order asymmetry parameter are shown in Table |^ [Q. An interesting feature of 
the results are that the wide resonance is consistently found to be lighter than the narrow 
resonance, and that the narrow plus skewed Gaussian model shows only a 20% deviation 
from the QCD value of the fifth order asymmetry parameter. 
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m (GeV) 


M (GeV) 


Width (GeV) 






square 


1.33 ±0.18 


1.23 ±0.18 


0.95 ±0.12 


0.6 ±0.13 


-0.11 GeV"^ 


gauss 


1.41 ±0.19 


1.23 ±0.15 


0.52 ±0.06 


0.49 ±0.13 


-0.13GeVi° 


skewed 


1.38 ±0.13 


1.06 ±0.21 


0.69 ±0.07 


0.44 ±0.04 


-0.19 GeV^" 


QCD 










-0.24 GeV^o 



Table 3: Resonance parameters obtained in the various two-resonance scenarios. The 
label "square" denotes the narrow plus square pulse model, "gauss" refers to narrow plus 
Gaussian resonance model, and "skewed" indicates the narrow plus skewed Gaussian model. 
The mass M denotes the state associated with the quoted width, which corresponds to the 
equivalent Breit-Wigner width for the Gaussian models. The quoted resonance parameters 
include uncertainties introduced by the QCD input parameters, except for the fifth-order 
residual moment combination which is obtained from the central values. 

7 Comparison of Distributed Resonance Strength 
Models 

All the distributed resonance strength models considered in the previous section lead to the 
excellent agreement between the QCD and phenomenological sides of the normalized GSR 
shown in Figure |[ and clearly resolve the discrepancy evident in Figure |] corresponding 
to the single resonance model [0. However, a measure of the agreement between the 
theoretical and phenomenological curves in Figure || result in a which is an order of 
magnitude smaller in the two-resonance models [Q. The combination of this result with 
the Table |^ observation of at least an order of magnitude disagreement with the QCD 
value of the third-order asymmetry is compelling evidence in favour of a two-resonance 
scenario for distributed resonance strength, with an upper bound on the heavier mass of 
about 1.6 GeV obtained from Eq. and Table |. 

The various two-resonance models all have similar values of x^, which does not provide 
a useful criteria to distinguish between these models. However, the narrow plus skewed 
Gaussian resonance model leads to the best (~ 20%) agreement with the QCD value of 
the next-order moment. In this scenario, the mass predictions and the pattern of a lighter 
broad resonance and a heavier narrow resonance is consistent with the identification of 
gluonium content in the /o(1370) and /o(1500). 

8 Conclusions 

in this paper, the formulation of Gaussian sum-rules has been reviewed. The key qualitative 
feature of Gaussian sum-rules is their enhanced sensitivity to excited states in comparison 
with Laplace sum-rules. Thus any resonance strong enough to stand out from the QCD 
continuum should reveal itself in a GSR analysis. The significance of normalized GSRs has 
been emphasized since they provide information independent of the the finite-energy sum- 
rule constraint that arises from the evolution of GSRs through the diffusion equation [^. 
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Figure 5: Comparison of the theoretical prediction for Nq'^^ (s,t,so) with the dis- 
tributed resonance phenomenological models. The r values used for the three pairs of 
curves, from top to bottom in the figure, are respectively r = 2.0 GeV^, r = 3.0 GeV^, and 
r = 4.0 GeV^. The almost perfect overlap between QCD and phenomenology is particu- 
larly impressive because there are no free parameters corresponding to the normalization of 
the QCD prediction and phenomenological models. 

Methods for obtaining phenomenological predictions from GSRs have been reviewed 
and applied to the p meson to demonstrate the predictive power of GSRs [Q. These tech- 
niques are easily adapted to a variety of resonance models, and through combinations of 
GSR moments, provide a natural criterion for testing the accuracy of the phenomenological 
model in comparison with QCD Q. 

The important role of instanton effects in relation to the low-energy theorem has 
been illustrated for the GSRs of scalar gluonic currents. This analysis provides additional 
support for the interpretation of instanton effects first developed for Laplace sum-rules 
An identical, and natural, partitioning of instanton effects into a continuum and LET-like 
contribution resolves discrepancies between the LET-sensitive and LET-insensitive sum- 
rules in both the Laplace and Gaussian sum-rules. 

A detailed phenomenological analysis of the k = GSR is presented because it is is 
less sensitive to QCD parameter uncertainties than the k = —1 (LET-dependent) GSR. 
The analysis clearly rules out a single narrow resonance scenario, and further indicates 
that gluonium resonance strength is spread over a broad energy region. In particular, 
evidence exists for two resonances below ^ 1.6 GeV with a significant gluonium content, 
with the lighter resonance having a substantial width [^. Such a scenario is consistent 
gluonium content of /o(1370) and /o(1500). 
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